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1 Introduction 

The forward-backward asymmetry A b r b in the production of bottom quarks at electron- 
positron colliders displays presently a substantial discrepancy between experimental 
measurement and theoretical expectation PQ. This theoretical expectation is deter- 
mined using the electroweak parameters obtained from a global fit to a number of 
different electroweak precision observables, including this forward-backward asym- 
metry itself. In turn, A\ b exercises a strong pull on the global fit, in particular 
towards larger masses of the Higgs boson. In particular, among the electroweak 
precision observables, A\ b is more sensitive on the Higgs mass than most other 
quantities. 

The present theoretical description of A\ b includes the fully massive next-to- 
leading order (NLO) electroweak [2] and fully massive NLO QCD H] corrections 
as well as the leading terms from the next-to-next-to-leading order (NNLO) QCD 
corrections |H] (see also 011!), which were obtained based on the massless approx- 
imation plus leading logarithmic mass terms. Given the substantial discrepancy 
between the experimental result and the theoretical expectation and the high im- 
pact on the Higgs mass determination, a more precise theoretical understanding of 
A b j b is clearly desired. 

At a future linear collider jH] , precision determinations of electroweak parameters 
will again involve the forward-backward asymmetries. In this setting, the top quark 
asymmetry A\ b , is experimentally accessible, and of high interest in understanding 
the interplay of quark mass generation and electroweak symmetry breaking. For 
a precise theoretical description of this asymmetry, inclusion of mass corrections is 
clearly mandatory. 

The NNLO QCD corrections to Ai b for massive quarks Q involve three classes 
of contributions: (1) the tree level matrix elements for the decay of a vector boson 
into four partons, (at least) two of which being the heavy quark-antiquark pair; (2) 
the one-loop corrected matrix elements for the decay of a vector boson into a heavy 
quark-antiquark pair plus a gluon; (3) the two-loop corrections to the decay of a 
vector boson into a heavy quark-antiquark pair. While the former two contributions 
can be obtained |Hj along the lines of the calculations of three jet production involving 
heavy quarks [TUJ HH E2] , the latter remain to be calculated. 

It is the aim of this paper (and of two companion papers) to contribute to the 
NNLO QCD corrections to A^ b for massive quarks by computing the virtual two-loop 
QCD corrections to the form factors of a massive quark. These form factors describe 
the full structure of the (Z*,7*) — > QQ vertex function, involving the vector and 
axial vector couplings of the vector boson. In the present paper, we derive the two- 
loop corrections to the vector form factors, while two following papers will discuss 
the parity-violating form factors at two loops. The two-loop corrections to the vector 
form factors were considered previously only in [13J, where the contribution from 
closed fermion loops was calculated. 

This paper is organized as follows. 

In Section |21 we give the notations followed throughout all the paper, defining the 
form factors and their expansions in terms of the coupling constant. In Section |3] we 





Figure 1: Tree- level and one-loop diagrams, involved in the calculation of the 
heavy-quark vertex form factors. The curly line represents a gluon; the double 
straight lines, quarks of mass m. The external fermion lines are on the mass-shell: 
p\= p\ = m 2 - The wavy line on the l.h.s. carries momentum Q = p± + P2, with the 
metrical convention Q 2 < when Q is space-like. 



give the virtual contributions to the one- and two-loop form factors before the sub- 
traction of the UV divergences. In Section0]we discuss in detail the renormalization 
and in Section El we give the UV-renormalized form factors at the one- and two-loop 
level. These results are presented for the case of choosing the renormalization scale 
\i equal to the mass m of the heavy quark. In the Subsection 15.31 the logarithms of 
the ratio m/fi, that are present if /j, ^ m, are explicitly given. Finally, Section El 
deals with the analytical continuation of our formulas above the threshold. 

2 The QCD Form Factors 

We call V£ c (pi,p2) the QCD vertex amplitude, corresponding to the decay of a 
virtual photon of momentum Q = P1+P2 into a quark and an anti-quark, of momenta 
Pi, p^ and colors c\ and c 2 respectively. The two fermions are on the mass-shell, 
p\= p\= m 2 , where m denotes the mass of the heavy quark in the on-shell scheme. 
Let us define the following two vectors: 

Q>=rt+i$, *=$-$., (i) 

such that Q 2 = S, where S is the cm. energy squared, and the dimensionless 
variable 

m m z 
Within QCD the V^ C2 (pi,p2) can be expressed in terms of two dimensionless 
scalar form factors -Fj(s), i — 1, 2, depending only on the dimensionless variable s of 
Eq. (|21), as follows: 

V£ C2 (pi,p 2 ) = u Cl (p 1 )T^ C2 ip u p 2 )v C2 (p 2 ) , (3) 






Figure 2: The two-loop vertex diagrams involved in the calculation of the form 
factors at order 0{a 2 s ). The curly lines are gluons; the double straight lines, quarks 
of mass m; the single straight lines, massless quarks and the dashed lines ghosts. 
The external fermion lines are on the mass-shell: p\ = "p\ = m 2 . The wavy line on 



the l.h.s. carries momentum Q 
Q is space-like. 



Pi +P2, with the metrical convention Q < when 
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where u cl (pi),v C2 (p2) are the spinor wave functions of the quark and the anti-quark, 
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[7^, 7"], and where: 



VQ = eQ Q} (5) 

Qq being the charge of the heavy quark in units of the positron charge e > 0. 

The indices C\, c 2 , in Eq. (JU), are the color indices: C\, c 2 — 1, ■ • • , N c , where N c 
is the number of colors. 

The extraction of each form factor -Fj(s) from Eq. (@J) can be carried out by the 
following general projector operators P^ . 
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performing a trace over the spinor and the color indices. The constants g ,j = 1, 2, 
are properly chosen such that: 

Tr (P^(m,p 1 ,p 2 )T^p 1 ,p 2 )) = F t (s). (7) 

Let us observe that since we work in a D-dimensional space (to regularize the diver- 
gences arising in the computation) the trace over the spinor indices is consistently 
performed in D dimensions as well. We use the convention of keeping the trace of 
the unit matrix equal to four also in D dimensions. 
The explicit values of the constants are: 
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The form factors are given as an expansion in powers of as/C2ir), where as is 
defined as the standard MS coupling in QCD (with Nf massless and one massive 
quark) at the renormalization scale \x: 
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The superscripts "1/" and "2/" stand for one- and two-loop contributions and the 
first term 1 in Fi(e,s, fi 2 /m 2 ) is the tree-level approximation. The subscript "i?" 
stands for "renormalized" , meaning that F>J(e, s,fi 2 /m 2 ) and F>J(e, s,fi 2 /m 2 ) 
come from the sum of the contributions of the virtual diagrams of Figs. ^ and 121 and 
the relative counterterms, shown in Fig. El for the subtraction of the UV divergences. 
In the following the expressions of the unsubtracted as well as the UV-renormalized 
form factors will be given at the one- and two-loop level. 

3 Unsubtracted Contributions 

In this Section we consider the contribution of the diagrams shown in Figs.^and|2to 
the expansions in Eqs. (J12I13JI . As explained in ^3JEHEI|, the form factors coming 
from the computation of the trace operation introduced in the previous Section are 
expressed in terms of several hundreds of scalar integrals. It is possible to express all 
these integrals as a combination of only 17 independent scalar integrals, called the 
Master Integrals (Mis) of the problem, by means of the so-called Laporta algorithm 
[T7] . using integration-by-parts identities [18J, Lorentz invariance ^H] and general 
symmetry relations. This reduction algorithm is performed exactly in D = (4 — 2e) 
dimensions |20j. Once the expression in terms of the Mis is found, we expand the 
result in powers of e around e = (D = 4), using the expansions of the Mis given in 
[HUTU] (the Mis are evaluated with the differential equations method [2*T| |2~2"1 I2H] ) . 
The result will be therefore given as a Laurent expansion in e where both UV- and 
IR-poles are regularized with the same parameter e. 

The form factors that we are going to present in this Section are given for space- 
like Q (S = Q 2 < 0) and they are expressed in terms of 1-dimensional harmonic 
polylogarithms (HPLs) J2U 123 of the variable: 



V-s + 4 + v^i ^/s + Am 2 + ^S ~ ~ 



Cf and Ca are the Casimir operators for the fundamental and adjoint represen- 
tation of the color gauge group SU(N C ) respectively: 

iV 2 -l 
Cf = H^vT' C a = N c , (15) 

where N c is the number of colors. Tr is the normalization factor of the generators 
of the fundamental representation: 

tr (t a t b ) = T R S ab = -5 ab . (16) 

Because all the calculations are done in D dimensions, we have to take into 
account for each loop a factor 

cwgy^rra^g)', (it) 

where /x is the mass scale of dimensional regularization, that we choose equal to the 
renormalization scale. 



3.1 One-Loop Unsubtracted Form Factors 

At the one-loop level only the diagram shown in Fig. ^ (b) contributes to the form 
factors Fi(e,s, fi 2 /m 2 ) and F 2 (e, s,/x 2 /m 2 ). We will write: 

^fU)( £ ' s '^) = C ^ {^^ (*>*)> with z = 1,2, (18) 
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3.2 Two-Loop Unsubtracted Form Factors 

At the two-loop level, all the 9 Feynman diagrams of Fig. |21 contribute to the form 
factors Fi(e,s, fi 2 /m 2 ) and F 2 (e, s,/i 2 /m 2 ). The total contribution comes from the 
sum of the diagrams involving only the heavy quark, Fig. |2] (a)-(e) and (g)-(i), and 
the diagrams in which a light quark runs in the internal loop, Fig.|2](f). If we consider 
Nf light quarks, the latter contribution is simply the contribution of diagram (f), in 
which the quark in the internal loop is considered as massless, multiplied by Nf. 
We will write: 
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4 4 4 



1 + x) ^ v ' (1-x) 

16 \ / 4 
X #l,0;x + 44 



(1-x) 2 



(1-x) (1 + x) 2 (1 + x) 



tf(l,0,-l,0;x) 



-4 + 



24 48 
+ 



38 



(1-x) (1 + x) 4 (1 + x) 3 (1 + x) 2 

18 \ / 4 4 
.|. 7r -l T )iJ(l,0,0;s)+(2-- -4 



;i+x) 

840 



336 



+ 



688 188 
4 



1 + x) 4 (1 + x) 3 (1+x) 2 (1+x) 



(1-x) 2 (1-x) (1 + x) 5 

16 

'#(l,0,0,0;x) 



"I 4 +7lfV7T+^ + 7T4Vdb) W 1 - - 1 - ^ 



2- 



1-x 1+x 



H(l,l,0:x 



, j-, w, . 



+ C F 7W 



169 34C(2) 22C(2) 22((2) 8((3) 

216 9(1 -x) 9(1 + x) 
3C(3) 8C(3) 8C(2) f I 



3(l+x) ' 3 
16 16 



9 3(1 -r) 

1 



16 



3 3(1 -x) 3(1 + x) 



1 — x) (1+x 
#(-1,-1, 0;x 



H(-l:x) 



+ ; r—, r #(-l,0;x)-, , 

1 9 9(l-x) 9(l+x)y v ' ' ; V3 3(1 -.<■) 



3(1 +x) 
175 

"27(1 + x) 



#(0;x) 



3(1 + x) 

22 
9(1 + x) 

4 
3(1 + x) 



#(0,-l,0;x)- 
#(0,0; x)- 



3 3(1 -x) 
28 34 



9 9(1 -x) 
4 4 



3 3 1-x 



#(0,0, 0;x) 



+ C F T R 



"223 4C(2) | 392C(2) 784((2) | 458((2) 



216 3(1 -x) 3(1 + x) 4 3(1 + x) 3 3(1 + x) 
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68C(2) i 8C(2 ) 4 ^ 3 ) 4 ^ 3 ) | 4 ^ 3 ) 



3(1 + x) v ' 3(1 -x) 3(1 + x) 

196 196 4C(2)/ 1 I 



f-^- 7^ T + 



9(l+x) 2 9(l+x) 3 \(l-x) (1+x) 

■ 1 ) jr (- 1 ")-I( 1 -(i^-(i^y) H( " 1, " 1,0ix) 

/4 4 

# -l,0;x + -- 



3(1 -x) 3(l+x)/ v ' ' ' \3 3(1 -x 
4 \ . . /409 C(2) 24C(2) 



3(1 + x)J v ' ' ' ' V 54 (i-^) (1 + +) 5 
60C(2) + 44C(2) 6C(2) C(2) 308 



(1 + x) 4 (1 + x) 3 (1 + x) 2 (1 + x) 27(1 - x) 
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+ 777^ — ^ " ^77—^ )H(p;x) 



9(1 +x) 3 3(1 + x) 2 27(1 + re) 

+ (5"3(T^)-3(TT^))^ '- 1 ' 0;I) 
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9 3(1 -x) 9(1 + x) 4 9(1 + x) 3 9(1 + x) 

74 7 24 

H(0,0;x)- U- 0/1 r + " 
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60 44 6 7 \ 
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240C(2) 508C(2) 270((2) 8((2) 
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- (T^y ) H C. -L - 1 - <* *> + ( 13 + (T^) + (if^i 
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R 
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20/3 12 18 
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2 1 6 
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90C(2) 180C(2) 231C(2) 51((2) 



+ 
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(1 + x) 4 (1 + x) 3 2(1 + x) 2 2(1 + x) 

23 6 9 41 

+ 
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6(1 -x) (1 + x) 3 (1 + x) 2 6(1 + x) 
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2(1 -x) 2 2(1 -x) (1 + x) 4 (1 + x) 3 2(1 + xf 
+ ™) H( - 1Mx)+ (J«?L_ 11C(2) 3 C (2) 



2 1+x 



(1-x) 3 4(1 -x) 2 8(1 -x) 



, 258C(2) _ 744C(2) + 1515C(2) _ 621((2) + 307C(2) 



(1 + x) 5 (1 + x) 4 2(1 + x) 3 2(1 + x) 2 8(1 + x) 
| 13C(3) _ 324C(3) | 810C(3) _ 635C(3) | 285((3) 



4(1 -x) (1 + x) 5 (1 + x) 4 (1 + x) 3 2(1 + xf 
13C(3) 259 



9 27 89 .. „ „ , 

■h — - + — . ; )H{0;x) 



4(l+x) 18(l-x) (1+x) 3 2(l+x) 2 9(l+x) 

9C(2) 27C(2) 21C(2) 90C(2) 225((2) 

2(1 -x) 3 4(1 -x) 2 8(1 -x) (1 + x) 5 (1 + x) 4 
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174<(2) + 36«2) + 21«2) ^ 



1 + x) 3 (1 + x) 2 8(1 + a)/ v ' ' Vl 1 -^) 

1 12 24 21 9 

l-x) + (1 + x) 4 ~ (1 + x) 3 + (1 + x) 2 ~ (1 + x 
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+ 777: r ~TZ tt + 



2(1 -x) 3 4(1 -x) 2 8(1 -x) (1 + x) 5 (1 + x) 4 
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+ 



(1 + x) 3 (1 + x) 2 8(1 + x) 

7C(2) 21C(2) 5C(2) 69C(2) 345((2) 



4(1 -x) 3 8(1 -x) 2 16(1 -x) (1 + x) 5 2(1 + x) 



1 



135C(2) 30C(2) 5C(2) 3 59 

(1 + x) 3 (1 + x) 2 16(1 + x) (l-x) 2 12(1 -x) 

24 39 7 115 A rr /n n N 



(1 + x) 4 (1 + x) 3 2(1 + x) 2 12(1 + x) 

1 3 25 12 30 

+ T77 T + 



(1-x) 3 2(1 -x) 2 4(1 -x) (1 + x) 5 (1 + x) 
48 42 25 \ 
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(l-x) 3 4(1 -x) 2 8(1 -x) (1 + x) 5 (1 + x) 
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4 3 11 



'1 + x) 3 2(1 + x) 2 16(1 + x) 
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|i?(0, 1,0; x) 
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70 45 11 \ , 
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(1 + x) 3 (1 + x) 2 2(1 + x) 

3C(2) 336C(2) __ 840C(2) __ 660C(2) _ 150C(2) 



(l-x) (1 + x) 5 (1 + x) 4 (1 + x) 3 (1 + x) 2 

3C(2) 24 36 12 \ 

+ JT+xj + (1 + x) 3 " (1 + x) 2 + Jl + x)) { ' ; X) 
24 48 40 16 \ 

1 + x) 4 " (1 + x) 3 + (1 + x) 2 " TT+x) J ( ' ' ' X) 
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+ 
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(l-x) (1 + x) 5 (1 + x) 4 (1 + x) 3 (1 + x) 2 
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[1 + x) 
+ C F T R N f 



#(1,0, 0,0; a;) 



4C(2) + 4C(2) 



3(1 -x) 3(1 + x) V 3 ! 1 -^) 

49 49 



+ 
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9(1 -x) 9(1 + x) 



#(0,0; x) 



#(0;x) 



+ C F T t 
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3(1 -x) 3(1 + x) 

2C(2) 136C(2) 272C(2) 132((2) 10C(2) 

3(1 -x) (1 + x) 4 (1 + x) 3 (1 + x) 2 3(l+x) 

68 68 



+ 



+ 



4 4 

f — , N )#(-!, 0;x) 



+ 



3(1 +x) 2 3(1 + x) ' V 3(l-x)3(l + x 

3C(2) , 24C(2) 60C(2) , 42((2) 3((2) 



2(1 -x) (1 + x) 5 (1 + x) 4 (1 + x) 3 (1 + x) 2 
3C(2) 44 124 62 



2(1 + x) 9(1 -x) 3(1 + x) 3 (1 
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9(1 + a) 
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#(0;x) 
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60 42 



3(1 -x) 3(1 +x) 4 3(1 + x) 3 

3 24 
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+ 



2(1 -x) (1 + x) 5 
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(1 + x) 4 (1 + x) 3 (1 + x) 2 2(1 + x) 
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1+x) 
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| 256C(2) | 9C(2) 69C 2 (2) | 207C 2 (2) 327( 2 (2) 



(1+x) 2 (1 + x) 20(1 -x) 3 40(1 -x) 2 80(1 -x) 
171C 2 (2) 855C 2 (2) 3291C 2 (2) 1323C 2 (2) 327( 2 (2) 
(1 + x) 5 ~ 2(1 + x) 4 + 10(1 + x) 3 ~ 20(1 + x) 2 ~~ 80(l+x) 
10C(3) 7C(3) 168C(3) , 348C(3) 200C(3) , 23((3) 



l-x) 2 (l-x) (1 + x) 4 (1 + x) 3 (1 + x) 2 (1 + x 
44 + 44 + ( 9C(2) 6C(2) 180C(2) 



(1 + x) 2 (1 + x) \(l-x) 2 (l-x) (1 + x) 4 
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(1 + x) 3 (1 + x) 2 (1 + x)/ l ,X) \(l-x) 
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(23) 



4 Renormalization 



The subtraction of the one-loop sub divergences and the two-loop overall divergence 
in the two- loop graphs of Fig.|2]is performed in a hybrid scheme: we renormalize the 
heavy-quark wave function and mass in the on-shell (OS) renormalization scheme, 
while the coupling as and the gluon wave function are renormalized in the mod- 
ified minimal subtraction (MS) scheme. The counterterm diagrams to add to the 
unsubtracted form factors given in Eqs. ()19|20|22|23jl . are shown in Fig. El 

4.1 One-Loop Counterterms 

The renormalization of the one-loop form factors and the subtraction of the one-loop 
sub divergences from the two- loop graphs shown in Fig.|2]require the renormalization 
constants Z „Mg(e), ^3 Ms( e )' ^ m os( e , m,fi 2 /m 2 ) and Z 2) o5(e, /i 2 /m 2 ), at the one- 
loop level. Once Z ^(e), Z 3 ^(e) and Z 2i os(e, fi 2 /m 2 ) are defined, the expression 
of Zip(e, fi 2 /m 2 ) is obtained using the Slavnov- Taylor identities. 

The expressions for the coupling and gluon wave function renormalization con- 
stants (the latter one in the Feynman gauge) are in the MS-scheme: 



(24) 



Z%(e) = -2* C(e) - (-C A - -T R (N f + l)) , 
9 ,msV I 2vr ; 4e V 3 3 V ; ') ' 

Z %^ = % c ®h(l CA -* TR ( N ' +1) ) ■ (25) 
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(a) 



(b) 



(c) 





(d) 




(g) 



(h) 





(i) 



(J) 



Figure 3: Counterterm diagrams. For the renormalization of the two-loop form 
factors we use diagrams (a)-(i). Diagram (j) is employed in the renormalization of 
the one-loop form factors. 
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For what concerns the renormalization of the heavy quark mass and wave func- 
tion, we need the following constants at the one- loop level, defined in the OS-scheme: 



,(10 



/'• 



Sm ( g [ e, m, „ 



r{ll) 



P 



7 K > ^ 1 



m — 6 e — - 

2vr w U 2 



-Tc(e) (£, 



2n 



m' 



C F (3 - 2e) 
~2~e(l-2e) 
C F (3-2e) 
~2~e(l-2e) ' 



(26) 

(27) 



The constant ZJ F (e,/i 2 /m 2 ), needed for the renormalization of the QQ-gluon 
vertex, is obtained from a Slavnov- Taylor identity as follows: 



IF I 6 ' ™2 



z (i0 



77T 



^w 



,(i0 



Z2 '° 5 ^'m 2 



^ 



+ 2 Z 3^M5 (e 



-^C(e)i 

2vr w 2e 



/'• 



Ca + -^ C. 



m' 



(10 

i,M! 

(3 ~ 2e) 
(1 - 2e) 



(28) 
(29) 



The renormalization of the UV divergences of the one-loop form factors, Eqs. (J19I20J1 . 
is straightforward. It is sufficient to add the counterterm of Fig. El (j), defined as: 




def 



v(10 ( f 
■^2,051 fc ' 



/' Z 



m z 




(30) 



where the constant Z\ os (e, fi 2 /m 2 ) is given by Eq. (|2*7jl . to the diagram of Fig.[T](b). 
The corresponding U V-renormalized form factors are given in Section 03 Let us note 
that J~2 (e, s) in Eq. (|20J1 is IR and UV finite. The counterterm diagram defined in 
Eq. (|3~U|l is, in fact, proportional to 7 M and affects only T\ (e, s). 

For the subtraction of the one-loop sub divergences from the two-loop diagrams 
of Fig. |2l we have to define as well the counterterm diagrams shown in Fig.|3j(a)-(h). 

The first two involve the constant 5m g(e,m, fx 2 /m 2 ): 



p-aph (a) in Fig. (J3J): 




= <5mU(e, m, — -I a ! m 

m V m z 




(31) 



graph (b) in Fig. (J3J): 




SmPQ 



= Srr%J[e, m, —? ) x | m 

m V m z 



,(10/ 




(32) 



The one-loop diagram multiplying 5m g(e,m, n 2 /m 2 )/m in Eq. ()31j) is defined 
as follows: 




mC F — 6(e 
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2 \ e 



m z 
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x / T) D k 



W 



[(pi + k) 2 — m 2 ] \{jp2 ~ k) 2 — m 2 ]k 2 



where the measure T) D k is such that: 

1 



J) u k 
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C(e) \fi 2 J J (2tt) 2 ( 1 ^) 



d D k 



and where: 



U^ = VQ7 a [^i+^+m][^i+^+m]7 M [/(f-^ 2 +m]7 a . 
The corresponding form factors are: 
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The one-loop diagram appearing in Eq. (J32)) has exactly the same form factors 
given in Eqs. (J37J38J) . Therefore, in the following we will refer to both diagrams 
using the picture of the first one only. 

The counterterm diagrams which involve Z2 os(e> ^ 2 / m2 ), Zip(e, n 2 /m 2 ), and 

loop vertex diagram: 



j 3MS (e) are defined as the product of the renormalization constants times the one- 



graph (c) in Fig. (jSJ): 




graph (d) in Fig. (JHJ): 




graph (e) in Fig. (jSJ): 



def (ii) / //■ . 
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def (10 / ^_ . 

^2,05 ^ C > m 2 J X 




(39) 
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(ii) d ±f _ z (iO / 



3,MS V 



graph (f) in Fig. Q: 




graph (g) in Fig. (jSJ: 




graph (h) in Fig. (JBJ): 




= Z -(^ 1 X 



" Z -(^l x 



(41) 




(42) 




(43) 




def 7 {Xl) ( ^_ . 



The one-loop vertex diagram, appearing in Eqs. 




(44) 



is defined as: 




2ir \m 2 J J [(pi + k) 2 — m 2 ][(p 2 — k) 2 — m 2 ]k 2 



where 



V M = VQ^l^+ft+mjYlft-fa + m}^ 



(45) 



(46) 



and the corresponding form factors are given in Eqs. (jU 



4.2 Two-Loop Counterterm 

In order to complete the UV-renormalization of the form factor F\ , we have also to 
subtract its value at s = 0, or, which is the same, to add the counterterm diagram 
shown in Fig. 01 (i). We need, therefore, the constant Z 2t os( e , fi 2 /m 2 ) at the two-loop 
level, that was computed in J2E1 121] • We use the result of [27| and we express it in 
terms of the renormalized MS coupling as, finding: 
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(47) 



The counterterm diagram is defined, therefore, as follows: 





'(20 ^ 2 



5 Renormalized Form Factors 

We report now the analytic expression of the UV-renormalized form factors at the 
one- and two-loop level, F\ R (e, s,// 2 /m 2 ) and F ii? (e, s,/i 2 /m 2 ), in the space-like 
region s = g 2 < 0, in terms of HPLs of the variable x already introduced in Eq. (|14jl. 



V-s + 4 



V-s + 4 + v^ ' 



The one-loop form factors are given up to the first term in the expansion in e, while 
the two- loop ones up to the finite part. 

The renormalization of the UV divergences is carried out in the hybrid scheme 
explained in the previous Section. Note that adding the diagrams of Figs. Q and 121 
and the corresponding counterterms we will have a non-trivial dependence on the 
renormalization scale /i, due to the fact that the virtual contributions coming from 
the one- and two-loop diagrams reported in Section El depend on the ratio (/z 2 /m 2 ) 2e , 
while in the counterterms we have a dependence on (fi 2 /m 2 ) 2€ (those calculated in 
the OS-scheme) as well as on (/i 2 /m 2 y (for the counterterms in the MS-scheme). The 
expansion of these ratios in powers of e generates terms proportional to log (/x 2 /m 2 ) 
and log (/j 2 / m 2 ). In this Section we give the one- and two-loop UV-renormalized 
form factors for // = m. For the terms proportional to log (fj, 2 /m 2 ) and log (/x 2 /m 2 ) 
we refer the reader to Section I5~3l 

5.1 One-Loop UV-Renormalized Form Factors 

The one-loop UV-renormalized form factors are recovered adding the diagram of 
Fig. ^ (b) and the corresponding counterterm: 

CTM(e,8)= Z$ 5 (e) x^ (49) 



We put /i = m and we define: 

^S?(e, s) = C(e) jg?(e, s) , with . = 1,2, (50) 
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5.2 Two-Loop UV-Renormalized Form Factors 

and (|48jl . we find, diagram- 



Adding the contributions of Eqs. (JH1JH2J) and 
matically, the following counterterm: 
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The two-loop UV-renormalized form factors are recovered adding the diagrams 
in Fig. Eland the counterterm in Eq. ()5Hj) . 
We put \x = m and we define: 
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5.3 Form Factors for /x ^ m 

In this Section we report the expressions for the renormalized form factors in the 
case we keep \i ^ m. 

At the one-loop level we do not have an explicit dependence on the logarithm of 
the ratio of the renormalization scale and the mass of the heavy quark, because an 
overall factor (/i 2 /m 2 ) e can be taken out: 
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where the functions T\ R (e, s) are given in Eqs. 

At the two-loop level, such a dependence results from the coupling constant 
renormalization, first appearing at this level. Factoring an overall (fi 2 /m 2 ) 2e , we 
have: 
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where the functions jF ii? (e, s) are given in Eqs. (|55|5fij) and the functions M\ (e, s) 

and Ay (s) can be derived from the renormalization group equation. 
Introducing (3q, the first coefficient of the QCD /3-function: 
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we can write 
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From these, the coefficients appearing in Eq. (|58jl can be readily read off: 
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6 Analytical Continuation above Threshold 



Eqs. (I51I52I55I56|) are written in terms of the "space-like" variable x, defined in 
Eq. (|IH) , for negative cm. energy squared S < 0. In particular, for S < 0, x is real 
and positive and varies from x = 1, when S = 0, to x = when S = — oo. 

The physical form factors, defined in the time-like region S = Q 2 > (in partic- 
ular above the physical threshold S > Am 2 , where an imaginary part appears) can 
be recovered by analytical continuation with the usual ie-prescription, i.e. giving a 
small positive imaginary part to S: S + ie. 

In so doing, if S > 0, but still S < Am 2 , the variable x becomes a phase factor: 
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where: 
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Above threshold, S > Am 2 , we define: 
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y/S-y/S- Am 2 



(6£ 



with y = 1 at S — Am 2 and y = at S = oo, and the continuation in x is performed 
by the replacement: 

x — ► —y + ie . (69) 

The real and imaginary parts of the form factors are defined through the rela- 
tions: 
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where s = S/m 2 . 

In the following two Sections we will give real and imaginary parts of the one- 
and two-loop analytically continued form factors for fi = m. The renormalization 
scale dependence follows from the pattern outlined in Section 15.31 



6.1 One-Loop Form Factors above Threshold 

As in Section El we write: 

F%\e, s) = C(e) *g>(c, a) with i = 1, 2 . (72) 

We have: 
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6.2 Two-Loop Form Factors above Threshold 

We define: 

F$(6,s)=C 2 (6)FW(6,s) With i= 1,2. 

We have: 
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The two-loop contribution to the heavy quark vector form factors arising from 
closed massive and massless fermion loops were computed previously in |13 i . In this 
work, the result for the massless fermion loops was obtained from a more general 
(unintegrated) result with different masses for the external and virtual quarks, by 
expanding in the mass of the virtual quark. 

Comparing to |13j . we fully agree on the form factors for the massive fermion 
loop, and on the leading logarithmic terms for the massless fermion loop. Owing 
to the different regularization procedures used in the latter case (small quark mass 
versus dimensional regularization), the finite terms differ. 



6.3 Threshold Expansions 

In this Section we provide the expansions of our results in the threshold limit S 
4m 2 (y — *■ 1 in the transformed variable). We define: 
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as the small parameter in which we expand. Keeping terms up to the zeroth order 
in /3, we have: 
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The above results are in perfect agreement with the results already known in 
the literature. In particular, the threshold limits of the 1-loop form factors are 
in agreement with Eqs. (9,10) of [28J (if we put C F = 1 in our results), once the 
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following replacement is carried out: 
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Consequently, we found also agreement with the 1-loop correction to the cross- 
section of e + e _ — > f f given in Eqs. (6,11,12) of [2%] . 

At 2-loop level, Eqs. (26,28) of [28J, concerning the abelian contributions match- 
ing our Cp terms, are written regularizing the IR divergences with a small mass 
for the photon. If we replace In (A/M) using Eq. (|92|). we are able to match the 
poles in 1/e with our expressions calculated directly in dimensional regularization. 
Nevertheless, we can not find agreement for the finite parts because of the differ- 
ences between the two regularization schemes. Eqs. (27,29) of |2E], instead, are in 
complete agreement with our CfTr terms. Let us note, however, that the contri- 
butions at two loops to the cross section e + e~ — > ff are IR finite and they can 
not depend on the regularization scheme. This is actually the case and we found 
complete agreement with Eqs. (31,32) of |25] . 

The 2-loop corrections in the threshold limit to the cross-section e + e _ — > QQ 
were calculated in |3TJ E2] : 
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where a^ is the tree-level cross-section and A^ 1 ) and A^ can be expressed, up to 
0(/3 2 ), in terms of the form factors as: 
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We found complete agreement with the results presented in these papers. 



6.4 Asymptotic Expansions 

In this Section we provide the expansions of our results in the limit S ^> m 2 (y — > 
in the transformed variable). Putting L = In (S/m 2 ) and keeping terms up to the 
second order in (m 2 /S), we have: 
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All the results of this Section can be obtained in an electronic form by down- 
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7 Summary and Outlook 

In this paper, we calculated the two-loop QCD corrections to the vector vertex form 
factors for heavy quarks. The result for the electric and magnetic form factors was 
obtained keeping the full dependence of the form factors on the mass of the heavy 
quarks, as well as on the momentum transfer, which is maintained arbitrary. 

The extraction of the form factors from the Feynman diagrams involved in the 
calculation was carried out by means of standard projector operators. Each form 
factor is expressed, in this way, as a combination of several hundreds of scalar 
integrals, whose UV and IR divergences are regularized, within the Dimensional 
Regularization procedure, by the same parameter D, dimension of the space-time. 
Using the Laporta algorithm, it was possible to reduce the problem of the calculation 
of all these integrals to the calculation of 17 master integrals, already present in the 
literature. 

The renormalization of the UV divergences was carried out in a hybrid scheme 
in which the coupling constant and the gluon wave function are renormalized in the 
MS scheme, while the mass and wave function of the heavy quark are renormalized 
in the on- shell scheme. 

The expressions of the unsubtracted as well as the UV-renormalized form factors 
are given in a closed analytic form as a Laurent expansion in e = (4 — D)/2. The 
coefficients of this expansion have a suitable representation in terms of 1-dimensional 
harmonic polylogarithms. The presence of poles in e in our results is related to the 
fact that IR divergences are still present. These divergences have to be canceled 
against the divergences arising from the real radiation, which in this paper was not 
taken into account. 

Besides being part of the full NNLO corrections to the forward-backward asym- 
metry of heavy quarks, the results presented in this paper can, on their own, already 
be used in a number of applications. 

An immediate point of interest is the behaviour of the inclusive heavy quark 
production cross section above the threshold. In the continuum the 47r-integrated 
cross section was computed to order a 2 s J2H] and, in view of top quark pair production 
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at a future linear collider, very detailed NNLO studies have been carried out in the 
threshold region (see 30J for a review), where the cross section is most sensitive on 
the top quark mass. The form factors derived here can be extrapolated at the two- 
loop level to their threshold values. We find complete agreement with all threshold 
results available in the literature [2E1 EU E2] • 

All available calculations of NNLO QCD corrections to inclusive heavy quark 
production at asymptotically large energies were made using the optical theorem 
(see [33] for a review), which avoids the explicit calculation of the two-loop form 
factors. Computing more differential quantities such as rapidity distributions re- 
quires the explicit knowledge of the two-loop form factors as well as of the massive 
(or leading-mass) single and double real radiation corrections. In the massless case, 
such calculations could be performed only very recently [31] , and a first step towards 
a fully massive calculation would certainly only consider the leading mass terms. In 
view of this application, we also provided the leading mass expansions of our full 
two-loop corrections to the form factors. 

The singularity structure of the two-loop heavy quark form factors could also 
provide insight into the generic singularity structure of two-loop integrals involving 
massive quarks. The corresponding structure of massless two-loop QCD amplitudes 
has been predicted from non-abelian exponentiation in [3J3 EH], an d proven very 
valuable in the calculation of massless two-loop four-point amplitudes [37] ■ For QCD 
amplitudes involving massive quarks, the singularity structure is only understood 
at the one-loop level at present [SB]- At the two- loop level, we observe that the 
infrared divergent contributions to the form factors which are proportional to the 
colour factor C\ exponentiate naively [3JJ], as already seen in the QED calculation 
|15j . The other colour factors contain explicit 1/e singularities which can not be 
explained from naive abelian exponentiation and deserve further investigation. 

Finally, this calculation demonstrates that it is possible to analytically compute 
two-loop vertex functions with at least one internal mass scale. These functions 
appear in a variety of applications, ranging from electroweak corrections to heavy 
quark physics, where first genuine two-loop results were obtained only very recently 
|4U| |4*T1 W2\ 14*3] , using the same methods (reduction to master integrals, differential 
equations, and basis of harmonic polylogarithms) as employed here. Many more 
results in this domain are yet outstanding, and the methods here can be clearly 
instrumental for their derivation. 

To compute the NNLO QCD corrections to the forward-backward asymmetry 
of heavy quarks, one also needs the two-loop corrections to the axial vector form 
factors. These corrections subdivide into two independent classes: anomalous and 
non-anomalous diagrams. The calculation of these is currently in progress, and 
results will be reported in two future publications. 
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